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ABSTRACT 

The strong shocks in young supernova remnants (SNRs) should accelerate cosmic rays (CRs) 
and no doubt exists that relativistic electrons are produced in SNRs. However, direct and con- 
vincing evidence that SNRs produce CR nuclei has not yet been obtained and may, in fact, be 
long in coming if current 7-ray observatories do not see an unambiguous pion-decay feature. Nev- 
ertheless, the lack of an observed pion-decay feature does not necessarily mean that CR ions are 
not abundantly produced since ions do not radiate efficiently. If CR ions are produced efficiently 
by diffusive shock acceleration (DSA) , their presence will modify the hydrodynamics of the SNR 
and produce morphological effects which can be clearly seen in radiation produced by electrons. 
We describe in some detail our CR-hydro model, which couples DSA with the remnant hydrody- 
namics, and the synchrotron emission expected for two distinct parameter sets representing type 
la and type 11 supernovae. Several morphological features emerge in radial profiles, including 
the forward shock precursor, which are observable with current X-ray observatories and which 
should definitively show if young SNRs produce CR ions efficiently or not. For the specific case 
of SN f 006, we conclude, as have others, that the extremely short X-ray scale heights observed 
near the outer shock argue convincingly for the efficient production of CR ions. 

Subject headings: Supernova Remnants, cosmic rays, shock acceleration, X-ray and radio emission 



1. INTRODUCTION 

Collisionlcss shocks in supernova remnants 
(SNRs) are believed to produce the majority of 
Galactic cosmic rays (CRs), at least up to the so- 
called "knee" near fO 15 eV (see Hillas 2005, for a 
recent review) . While there is little doubt from the 
synchrotron interpretation of radio observations 
that young SNRs produce GeV electrons, and this 
is probably true for TeV electrons as well from 
the interpretation of nonthermal X-rays, there is 



as yet no unambiguous direct evidence that SNRs 
produce relativistic ions. This is somewhat para- 
doxical considering that the observed electron to 
proton ratio in CRs is ~ O.Of and virtually all 
models of diffusive particle acceleration in col- 
lisionless shocks, the most cited mechanism for 
producing CRs, predict that ions receive far more 
energy than electrons (see, for example, Baring 
et al. f999, and references therein). Relativistic 
electrons, of course, radiate far more efficiently 
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than do ions, leaving open the possibility that a 
large majority of the energy in relativistic particles 
in SNRs lies in hard to see ions. In this paper, 
we model SNR evolution coupled with the effi- 
cient production of CRs (our so-called CR-hydro 
model, e.g., Ellison et al. 2004) and make a num- 
ber of predictions for the synchrotron emission 
from electrons which will be influenced by the 
presence of otherwise unseen relativistic ions. For 
a recent summary of observations and models of 
synchrotron emission in SNRs, see Cassam-Chena'i 
et al. (2005) which addresses many of the issues 
discussed here using a self-similar approach. 

In order to power CRs, the shocks in SNRs 
must be capable of placing ~ 10% of the super- 
nova (SN) explosion energy into relativistic ions 
over their lifetime (e.g., Drury 1983; Blandford & 
Eichler 1987). In fact, the strong shocks in young 
SNRs may be far more efficient than this (e.g., El- 
lison 2000; Hughes et al. 2000; Decourchelle et al. 
2000) and place enough energy in relativistic par- 
ticles so that nonlinear feedback effects modify 
the shock structure, the evolution of the remnant, 
and the radiative properties (e.g., Berezhko et al. 
1996; Decourchelle et al. 2000). As we show be- 
low, structural changes produced by DSA trans- 
late into changes in synchrotron emission that are 
large enough to be investigated with modern, high 
spatial resolution, radio and X-ray observatories. 
In particular, we calculate the synchrotron emis- 
sion profiles for typical shell-type la and II super- 
nova parameters and show how these profiles pro- 
vide important constraints on the underlying par- 
ticle acceleration mechanism and magnetic field 
structure. 

Particle acceleration influences the SNR evo- 
lution because relativistic particles produce less 
pressure for a given energy density than do non- 
relativistic particles. 1 Therefore, when relativistic 
particles are produced and/or energetic particles 
escape from the shock system, the shocked gas be- 
comes more compressible, i.e., it acts as if it has a 
softer equation of state and the remnant hydrody- 
namics are modified. The softer effective equation 
of state means that compression ratios well in ex- 
cess of four can be produced in non-radiative, col- 

x This follows since the ratio of specific heats, T, decreases as 
particles become relativistic and the pressure P = (T — l)e, 
where e is the energy density. 



lisionless shocks (e.g., Blandford & Eichler 1987; 
Jones & Ellison 1991), and since the energy going 
into relativistic particles is drawn from the shock- 
heated thermal population, the temperature of the 
shocked gas can be much less than that expected 
from test-particle shock acceleration (e.g., Ellison 
2000; Decourchelle et al. 2000). In addition to 
modifying the evolution and the temperature of 
the shocked gas, changes in the compression of the 
fluid should result in changes in the compression of 
the magnetic field implying that synchrotron emis- 
sion from relativistic electrons will vary strongly 
with the efficiency of DSA and the orientation and 
strength of the magnetic field. 

Perhaps the most important morphological as- 
pect of this CR-hydro coupling is that the ratio 
of the forward shock radius, -Rfs> to the radius 
of the contact discontinuity, Rqd, may be much 
less than in the test-particle case (see Decourchelle 
et al. 2000; Ellison et al. 2004). If, as is generally 
believed, shocks put far more energy into accel- 
erated ions than electrons, it is the efficient pro- 
duction of cosmic ray ions that reduces Rfs I Rcd 
from test-particle values. However, since the inter- 
action region between the forward shock (FS) and 
the contact discontinuity (CD) can be sometimes 
estimated or determined with modern X-ray tele- 
scopes (SN 1006 is an example where the CD is not 
seen), radiating electrons can reveal the presence 
of these otherwise unseen relativistic ions. 

Another clear morphological prediction from ef- 
ficient DSA discussed below is that radial pro- 
files of X-ray emission will be strongly peaked and 
form sheet- like structures at the FS. This effect 
comes largely from the large shock compression ra- 
tios which compress the magnetic field behind the 
FS and result in severe radiative losses for elec- 
trons producing X-rays. Without efficient particle 
acceleration, the radial profiles of X-rays will be 
smoother and more closely resemble those for ra- 
dio emission. 

In addition to the radio and X-ray profiles in 
the interaction region between the CD and the 
FS, we calculate the emission in the FS precur- 
sor. We show that the structure of the X-ray pre- 
cursor depends strongly on assumptions made for 
the magnetic field compressibility. If the magnetic 
field is compressed substantially at the FS, as is 
likely, the ratio of X-ray intensity immediately up- 
stream of the shock to that at the FS drops dra- 
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matically. In this case, line-of-sight projection ef- 
fects produce profiles that are fully consistent with 
the extremely short scale heights seen in SN 1006 
by Bamba et al. (2003) or Long et al. (2003), even 
though TeV electrons with long diffusion lengths 
are present. We conclude for this particular rem- 
nant, as did Berezhko et al. (2003), that CR ions 
are being efficiently produced and their presence is 
revealed by radiating electrons. We note that the 
strong magnetic fields we describe at the FS are 
produced by compression not from magnetic field 
amplification resulting from cosmic-ray streaming 
instabilities, such as predicted by Bell & Lucek 
(2001). Magnetic field amplification at the FS is 
not included in our model. 

2. CR-HYDRO MODEL 

We calculate the hydrodynamic evolution of a 
SNR coupled to efficient DSA with a radially sym- 
metric model described in detail in Ellison et al. 
(2004) and references therein. We do not consider 
CR production at the reverse shock since we as- 
sume the magnetic field in the ejecta is the frozen- 
in field from the SN progenitor and, as such, will 
be too small to produce significant particle accel- 
eration or non-thermal emission without large en- 
hancement factors (see Ellison et al. 2005, for a 
discussion of efficient DSA at reverse SNR shocks). 

Any realistic model of a SNR will have sev- 
eral parameters for both the environment and the 
physical processes controlling the evolution and 
particle acceleration. Here, we concentrate on 
changes in the SNR evolution and emission pro- 
duced by CR production, and choose two fairly 
distinct models as prototypes, one with parame- 
ters typical of type la SNe and the other with pa- 
rameters likely those of type II SNe. These models 
differ by the initial density profile in the ejecta 2 
and the density and magnetic field profiles in the 
ambient medium. Within these models, we inves- 
tigate the effects of varying the CR production 
efficiency and the magnetic field structure. 

2 Since we don't consider acceleration at the reverse shock, 
the different ejecta composition in type la and type II SNe 
is not important. For a discussion of how composition 
might influence DSA, sec Ellison ct al. (2005). 



2.1. Type la prototype 

For our type la prototype, we assume the den- 
sity profile of the ejecta material is exponential 
(Dwarkadas & Chevalier 1998), the total ejecta 
mass is M c j = 1.4 M©, the explosion energy is 
Esn = 10 51 erg, and a uniform ambient medium 
density, n p , with a temperature of T = 10 4 K. 
Here, n p , is the proton number density and we as- 
sume there is an additional 10% contribution of 
helium nuclei. We assume the magnetic field in 
the interstellar medium (ISM), B^, is also con- 
stant and take B = 10~ 5 G as a default value. 
We typically view our type la models at an age 
isNR = 400 yr, similar to the age of Tycho's SNR, 
when the shock speed is roughly 4000 km s _1 . 

2.2. Type II prototype 

For our type II prototype, we assume the initial 
density profile of the ejecta material is a power 
law in radius, p c j oc r~ n , with a constant density 
plateau region at small radii (e.g., Arnett 1988). 
We take n = 9 in all of our type II models. For 
the total ejecta mass we take M e j = 2M Q , and 
the explosion energy is set to i^sN — 3 x 10 51 erg 
(Laming & Hwang 2003; Chevalier & Oishi 2003). 
The density of the pre-SN wind is taken as p w = 
Ar~ 2 , where A = dM/dt/(4 tt v w ), dM/dt is the 
mass loss rate, and v w is the wind speed (both 
assumed constant). We use typical values w w = 20 
km s" 1 , dM/dt = 2 x 10~ 5 M Q yr" 1 (Chevalier & 
Oishi 2003), and take a constant wind temperature 
T w = 10 4 K. 

Following Chevalier & Luo (1994), we assume 
the unshocked magnetic field in the pre-SN wind 
is 

B (r) = (a w v w dM/dt) 1/2 /r , (1) 

or 



Bo(r) 




where cr w is the constant ratio of magnetic field en- 
ergy density to kinetic energy density in the wind. 
This expression assumes that the magnetic field 
is frozen in the constant stellar wind and is only 
valid in the equatorial plane for distances r, much 
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greater than the radius of the pre-SN star. 3 Off 
the plane, B(r) will fall off more rapidly than 1/r, 
but we ignore this effect in our spherically sym- 
metric models. The value of a w for stars other 
than the sun is not well known but, for concrete- 
ness, we take a w = 0.1. We typically view our type 
II models at isNR = 400 yr, to match our type la 
models and for comparison to SNR Cassiopeia A 
(Cas A), when the shock speed is roughly 6000km 

2.3. Acceleration model 

For the diffusive shock acceleration process, 
we use the algebraic model of Berezhko & Elli- 
son (1999) and Ellison et al. (2000) where the 
injection efficiency is parameterized and the su- 
perthermal spectrum, f(p), is a broken power law, 
/pl(p) 5 with an exponential turnover at high mo- 
menta, f(p) oc /pL(p)exp(-p/p max ). The alge- 
braic model solves the nonlinear DSA problem at 
each time step of the hydro simulation given the 
shock speed, shock radius, ambient density and 
temperature, and ambient magnetic field deter- 
mined in the simulation. With the accelerated 
distribution, an effective ratio of specific heats is 
calculated and used in the hydrodynamic equa- 
tions, completing the coupling between the two 
(see Ellison et al. 2004, for a full discussion). 

The injection parameter, rj^, is the fraction of 
total protons injected into the DSA process and 
values ?7inj > 10~ 4 typically yield efficient particle 
acceleration rates where 10% to 99% of the avail- 
able energy flux goes into rclativistic protons. 4 
The maximum momentum, p max , is determined 
by setting the acceleration time equal to the SNR 
age tsNR or, by setting the diffusion length of 
the highest energy particles equal to some frac- 

3 Equation (1) only applies if the forward shock has not 
reached the stellar wind termination shock. We assume the 
forward shock is within the bubble in all of the examples 
discussed here. 

4 Note the difference between the fraction of protons injected 
into the acceleration process, ?7inj, and the acceleration ef- 
ficiency. The acceleration efficiency is the fraction of total 
energy flux going into relativistic particles including all ions 
and electrons. Given »7i n j and the other shock parameters, 
the electron spectrum is determined with two additional 
parameters, the electron to proton ratio at relativistic en- 
ergies, (e/p) re i, and the electron to proton temperature ra- 
tio immediately behind the shock, T e /T p (see Ellison et al. 
2000, for a full discussion). 



tion, / s k, of the shock radius i? s k, whichever gives 
the lowest p max (see, for example, Baring et al. 
1999). In all of the models presented here we take 
fsk = 0.05. We assume Bohm diffusion so that 
the scattering mean free path, A, is on the or- 
der of the gyroradius, r g , i.e., A = r] m i p r g with 
??mfp = 1 and r g = pc/(qB). Here, p and q are 
the particle momentum and charge, respectively, 
B is the magnetic field at the acceleration site, and 
c is the speed of light. Note that while our esti- 
mate of Pmax requires a specific assumption for the 
mean free path, the acceleration model itself only 
assumes that the scattering mean free path is a 
strongly increasing function of momentum. In the 
absence of radiative losses, the maximum kinetic 
energy particles receive in DSA depends only on 
the particle charge and p max is the same for pro- 
tons and electrons as long as both are relativistic. 

2.4. Synchrotron emission and losses 

As the forward shock overtakes fresh ambient 
medium material, the shock accelerates these par- 
ticles and produces a nonthermal distribution as 
described in detail in Ellison et al. (2004) and El- 
lison et al. (2005). 5 Once the particle distribution 
is produced in a shell of material at the shock, it is 
assumed to remain in that shell as the shell con- 
verts and evolves behind the shock. During the 
evolution, particles experience adiabatic and syn- 
chrotron losses and these losses are calculated as 
in Reynolds (1998). 

In calculating the synchrotron emission and 
losses, we evolve the magnetic field as described, 
for example, in Reynolds & Chevalier (1981) or 
Reynolds (1998). Consider a fluid element which 
is now at position r with density p(r). At an ear- 
lier time, this fluid element was shocked at a posi- 
tion ri where the density immediately behind the 
shock was p2- The radial and tangential compo- 
nents of the field immediately behind the shock 
at ti, were B^ r and i?2t, respectively. If the mag- 
netic flux is frozen in the fluid, the field at the 
downstream position, r, is given by 




For the magnetic field configuration across the 



5 We ignore pre-existing CRs and inject and accelerate only 
thermal particles over taken by the shock. 
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shock, we assume either that B 2 = Bq, as in a 
parallel shock, or that the field is fully turbulent 
upstream and, following Volk et al. (2002), set the 
immediate downstream magnetic field 

B 2 = ^1/3 + 2^/3 B Q , (4) 

where r to t is the shock compression ratio. 6 Note 
that B 2 does not include any amplification effects 
such as described by Bell & Lucek (2001). Using 
B(r) obtained in cq. (3), the evolution of the elec- 
tron distribution under combined adiabatic and 
synchrotron losses is calculated and, at the end of 
the simulation, the synchrotron emission in each 
shell is determined as in Baring et al. (1999). 7 

In Fig. 1 we show electron momentum phase- 
space distribution functions, f(p), for a type la 
SNR model discussed more fully in Section 3 be- 
low. In each panel, the dashed curve is the dis- 
tribution calculated immediately after production 
at the age indicated (i.e., at i s hock) an d the solid 
curve is this distribution at the end of the simu- 
lation (i.e., at ^snr = 1000 yr) after experiencing 
adiabatic and synchrotron losses. In the top two 
panels, the dot-dashed curves show the electron 
distribution at tsNR = 1000 yr when only adia- 
batic losses are included. 

The shock accelerated distribution, before 
losses, is a broken power law above a thermal dis- 
tribution with an exponential cutoff at the max- 
imum momentum (i.e., eq. 12 in Ellison et al. 
2000, with a = 1). Adiabatic losses affect all 
particles (shifting the entire distribution to lower 
momenta, i.e., p oc p 1 ^ 3 ), while synchrotron influ- 
ence mainly the highest energy electrons. For the 
parameters of this model, the highest momentum 
electrons accelerated at early times are strongly 
depleted and a distinct synchrotron bump is ob- 
served just below the sharp maximum momentum 
cutoff. The heavy-weight dotted curve in the bot- 
tom panel of Fig. 1 is the electron distribution at 
the end of the simulation summed over the inter- 
action region between the contact discontinuity 
and the forward shock. For comparison we show 

'Here and elsewhere the subscript (2) indicates values im- 
mediately ahead of (behind) the shock. 

In calculating electrons losses, we include inverse-Compton 
losses off the cosmic microwave background radiation as 
described in Baring et al. (1999). For protons, radiative 
losses are unimportant for typical SNR magnetic fields. 



with the light-weight dotted curve the summed 
proton distribution at the end of the simulation. 
For this example, the electron to proton ratio at 
relativistic energies, (e/p) rc i, is set to 0.01, similar 
to that of Galactic cosmic rays, and the electron to 
proton temperature ratio immediately behind the 
shock, T e /T p , is set to 1 (see Ellison et al. 2000, 
for fuller discussion of these parameters) . The dif- 
ference between the electron and proton spectra 
in the bottom panel of Fig. 1 illustrates how DSA 
typically puts far more energy into protons than 
electrons. 

2.5. Upstream precursor 

The algebraic acceleration model of Bcrczhko 
& Ellison (1999) doesn't explicitly include the ge- 
ometry of the shock precursor. However, we can 
estimate the precursor upstream of the forward 
shock in the following way. At any particular 
time, the proton distribution in the outer most 
shell, f p (p), produces the precursor. We assume 
that the protons of momentum p in this distribu- 
tion "feel" a flow speed u(z) and magnetic field 
B(z), where z is the diffusion length, Ld(p), mea- 
sured upstream from the FS. The diffusion length 
Ld{p) = k{p)/u{z), where n = Aw/3 is the diffu- 
sion coefficient, v is the particle speed, and u(z) 
is the flow speed at z measured in a frame at rest 
with the shock. 

We use information from f p (p) to estimate u(z) 
and B(z) and obtain Ld{p). Because of shock 
smoothing, the compression ratio in the FS that 
produced f P (p) ranges from the subshock compres- 
sion, r su b, felt by protons with the superthermal 
injection momentum p; n j, to the overall compres- 
sion, r to t, felt by protons withp max . Intermediate 
values of compression, r(p), felt by protons or elec- 
trons with momentum p between p m j and p max , 
can be estimated with a linear extrapolation be- 
tween r(p) and log (pv), i.e., 

r(p) = r suh + G(p) (r to t - r suh ) , (5) 

where pv is proportional to the diffusion length 
and 

g(p) = .yf'f;^ . (6) 

log (pw) m ax ~ fog tfWjinj 

Here (pv) max = Pmax c, (fw)inj = Pinj «inj, and i; inj 
is the particle speed corresponding to pi n j. Note 
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that since Pi n j and r/mj combine to give a single 
free injection parameter, the specific value of pi n j 
is unimportant for the results discussed here (see 
Blasi et al. 2005, for recent work on injection in a 
semi- analytic, nonlinear DSA model). 

With equation (5), we estimate the flow speed 
felt by a particle with momentum p as 

u(z) = VJ^ , (7) 

and the magnetic field felt by this particle is either 
B{z) = B (8) 

or 

depending on if the magnetic field is compressed 
in the precursor (as in Eq. 4) or not. Here V sk is 
the forward shock speed in the rest frame of the 
SN. 

Given u{z) and B(z), the diffusion length of an 
electron can be determined and, in a fashion sim- 
ilar to Reynolds (1998), we assume that electrons 
of momentum p are distributed upstream from the 
shock such that 

/ e (p, z) = f e (p,0)exp{-z[l/L D (p) + l/(f sk R FS )]} , 

. (10) 

where / c (p, 0) is the electron distribution in the 
outer most shell (z = 0) at the end of the simu- 
lation and ./sk-RFS sets the distance ahead of the 
shock where particles freely leave the system. The 
electron distribution, / c (p, 0), contains the effects 
of synchrotron and inverse-Compton losses which 
occur during acceleration. 

The above relations are approximations in that 
they ignore the precise form for the smooth pre- 
cursor flow speed. However, we have verified 
that the precursor emission is relatively insensi- 
tive to this smoothing and that our approxima- 
tions adequately describe the spatial dependencies 
important for predicting the synchrotron precur- 
sor. Typical results are shown in Fig. 2 where the 
solid curves are for compressed B and the dotted 
curves are for uncompressed B. 



3. RESULTS 

3.1. Radial emission 

Using the parameters for our type la prototype, 
we plot in Fig. 3 the synchrotron emission as a 
function of radius for one radio (1-1.4 GHz; solid 
curves) and two X-ray bands (0.1-1 keV dashed 
curves; 1-10 keV dotted curves). We present four 
models, two with ?7i n j = 10~ 3 , which produces very 
efficient DSA with nearly 100% of the energy flux 
crossing the shock going into relativistic particles, 
and two with ?7i n j = 10 , which yields essentially 
a test-particle result with less than 1% of the en- 
ergy flux going into CRs and where the influence 
of shock accelerated protons on the hydrodynam- 
ics is small. For each ?7i n j we show a case with a 
compressed field (labelled B comp.) and one with 
uncompressed field either in the shock or the pre- 
cursor (labelled B2 = Bq). In the compressed field 
case, we assume, as in Bcrezhko ct al. (2002), that 
the magnetic field is fully turbulent upstream of 
the shock and is compressed in the precursor as 
described by equation (9). The curves are nor- 
malized to each energy band's flux at the forward 
shock. 8 

Fig. 4 shows similar results for our type II pro- 
totype where, as in Fig. 3, the emission is viewed 
at t SNR = 400 yr. 

Comparing these figures, we note the following: 
(i) The two SN types have very similar profiles 
at least for the parameters used here. One no- 
ticeable difference occurs for the r?; n j = 10~ 5 cases 
where the type II radio profiles are flatter than the 
type la profiles. Later, in association with Fig. 9, 
we show in more detail that changes in 7/i n j and 

The results of the CR-hydro model, at early times, depend 
on the initial conditions which, unavoidably, are somewhat 
arbitrary. The initial conditions, in turn, influence the 
emission at the CD seen in Figs. 3 and 4. For all of the 
results presented here, the simulation is started at a time 
to = 10 yr with the initial ejecta speed varying linearly with 
radius from zero to a maximum speed Vmax = 0.1c. The 
initial maximum radius of the ejecta is set by Vm ax and to 
and the early stages of the simulation, and therefore the 
synchrotron emission at the CD, depend on Vmax an d to- 
Of course, the later evolution of the SNR is nearly indepen- 
dent of the starting conditions, as long as the total kinetic 
energy and ejecta mass stay the same. Since the X-ray 
emission is dominated by losses at the CD, it is only the 
radio emission at the CD that depends strongly on Vm J ax 
and to- F° r a full discussion of the start up conditions for 
the CR-hydro model, see Ellison ct al. (2004). 
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other parameters influence the SN types rather dif- 
ferently and may offer help in distinguishing the 
types. 

(ii) In the interaction region between the con- 
tact discontinuity and the forward shock, the X- 
ray synchrotron falls off more rapidly than the ra- 
dio emission. As mentioned in discussing Fig. 1, 
the electrons producing the radio emission suffer 
only adiabatic losses, while the higher energy elec- 
trons producing the X-rays suffer adiabatic losses 
combined with synchrotron and inverse-Compton 
losses. In Fig. 5 we show profiles for the 1-10 
keV band with no losses (solid curve), with just 
adiabatic losses (dashed curve), and with adia- 
batic plus radiative losses (dotted curve). Since, 
for typical SNR parameters, the nonthermal X-ray 
emission comes from the exponential part of the 
electron spectrum, the X-ray emission will be ex- 
tremely sensitive to changes in the spectrum com- 
ing from any type of loss mechanism. 

(iii) The radio emission can have a secondary 
peak at the CD, while the X-ray emission, with 
synchrotron losses, always drops precipitously at 
the CD. As just mentioned, the radio emission at 
the CD is sensitive to the starting conditions of 
the hydro model but, in any case, the secondary 
peak is less noticeable in projection as we show 
below. 

(iv) With efficient DSA and a compressed mag- 
netic field (top panels of Figs. 3 and 4), the X-ray 
fall-off is extremely rapid and the X-ray emission 
can appear as an extremely thin sheet at the FS. 

(v) The precursor outside of the FS falls slowly 
if the magnetic field is not compressed at the 
shock, but drops sharply immediately upstream of 
the shock when B is compressed, with or without 
efficient DSA (top two panels in Figs. 3 and 4). 
The sharp drop due to the compressed field will 
make the X-ray precursor faint and difficult to de- 
tect compared to the emission at the FS. Without 
compression, the precursor should be observable, 
providing an important diagnostic for the mag- 
netic field structure. Note that the radio precursor 
has an extremely short upstream diffusion length 
for all cases and will not be detectable if the diffu- 
sive length scale is anywhere near as small as we 
predict. 

(vi) Comparing the ?7i n j = 10~ 3 panels against 
the ?7i n j = 10~ 5 panels in Figs. 3 and 4 shows 



that the distance between the CD and the FS is 
nearly a factor of two greater in the test-particle 
case than with efficient DSA. Since the limit of 
the shocked ejecta gives an idea of the position of 
the CD, Rfs/Rgd is measurable in several young 
SNRs with Chandra and XMM-Newton, making 
this morphological difference a powerful diagnos- 
tic for efficient DSA. 

In Fig. 6 we show the magnetic field structure, 
at isNR = 400 yr, in the transition region between 
the CD and FS for our two prototypes with com- 
pressed B and rji n j = 10~ 3 . The numbers at spe- 
cific points on the curves indicate the compression 
ratio, r to t, at the FS at the time that particu- 
lar parcel of gas was shocked. It is notable that 
ftot 3> 4 in all cases. The difference in r to t be- 
tween the two models comes about mainly from 
the lower magnetic field in the pre-SN wind for 
the type II model which results in larger compres- 
sion ratios. The end of the curves, marked with 
an open circle, show the immediate upstream, un- 
shocked magnetic field, -Bo, at £snr- For type la, 
Bq = 10 [iG and is independent of time, while for 
type II, Bo(r) falls off with radius as in equa- 
tion (I) and at £snr = 400yr is ~ 1.5 [iG. A 
thorough discussion of the influence magnetic field 
strength has on r to t is given in Ellison et al. (2005). 

3.2. Line-of-sight projections 

In Fig. 7 we show line-of-sight projections for 
some of the results shown in Fig. 3. Even in pro- 
jection, it is clear that the radio emission falls off 
less rapidly behind the FS than the X-ray emis- 
sion. Projection has little effect on the upstream 
precursor so the large differences seen in Fig. 3 
with and without magnetic field compression are 
similar in projection. The decrease in i?Fs/i?CD 
for efficient particle acceleration is less obvious in 
projection but, since the CD generally shows up 
via thermal X-ray emission, i?ps /Rcd remains an 
important diagnostic for the presence of efficient 
CR ion acceleration. Line-of-sight projections of 
the results shown in Fig. 4 are similar. 

An important feature that is in the line-of-sight 
projections and not in the radial profiles is the off- 
set of radio and X-ray peaks at the FS. In Fig. 8, 
the projections for the type la models of Fig. 3 
with compressed magnetic fields are plotted as a 
fraction of the FS radius. With or without ef- 
ficient DSA, the radio peak (solid curve) occurs 
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inside the X-ray peaks. Behavior such as this is 
observed in several SNRs including G347 (Lazen- 
dic et al. 2004), Kepler (DeLaney et al. 2002), Ty- 
cho (Decourchelle et al. 2001), and Cas A (Long 
et al. 2003). We note, however, that there is an- 
other radio peak coincident with the X-ray peak 
in Tycho (e.g., Dickel et al. 1991). For the efficient 
acceleration case (top panel) , the two X-ray peaks 
are also well separated. Note also that because 
of projection effects, the maximum emission oc- 
curs inside of the FS. As emphasized by Berczhko 
et al. (2003), care must be taken not to interpret 
the peak emission as the position of the FS, as 
done by Bamba et al. (2003) for SN 1006. The ac- 
tual upstream precursor is indicated in Fig. 8 with 
a "P." 

In Fig. 9 we compare the line-of-sight 1-10 keV 
X-ray projections for both type la and type II pro- 
totypes calculated with different DSA injection ef- 
ficiencies. While the absolute normalization is ar- 
bitrary, the curves show the correct relative nor- 
malization between the various models and, as ex- 
pected, the test-particle cases with Tfrnj = 10~ 5 
have lower absolute emissivities. In both pan- 
els, the solid curves have 7/inj = 10~ 3 , the dashed 
curves have 7/inj = 10~ 4 , the dotted curves have 
Vinj — 10~ 5 , and all models have magnetic field 
compression (note the different vertical scales in 
the two panels). For both SN types, the ra- 
tio i?FS /-Rcd increases noticeably as the accelera- 
tion becomes less efficient, but Rfs/Rcd increases 
somewhat more rapidly for type II SNRs. Also, for 
both SN types, the morphology of the X-ray emis- 
sion varies strongly with rji n y. for efficient DSA, 
there is a pronounced peak at the rim, while the 
emission is much broader for inefficient DSA. This 
difference offers another important diagnostic for 
efficient DSA. 

In Fig. 10 we keep all parameters of our ?7i n j = 
10~ 3 type II model constant except the wind 
speed, w w , and the mass loss rate, dM/dt. In 
the top panel, w w = 20 km s _1 and dM/dt varies, 
as indicated, and the light-weight dashed curve 
has ?7i n j = 10 -4 ; all other curves in Fig. 10 
have r/inj = 10~ 3 . As dM/dt increases, there is 
an increase in Rfs/Rcd indicating, among other 
things, that self-similarity is no longer a good ap- 
proximation at tsNR = 400 yr. In the bottom 
panel, dM/dt = 2x 10~ 5 M yr -1 and v w is var- 
ied as indicated. Now, the profiles are relatively 



insensitive to the changes in w w , suggesting that 
self-similarity does apply. 

In considering Figs. 9 and 10 it's important to 
note that while Rfs/Rcd is reduced substantially 
with efficient CR production in type la SNRs, 
values of Rfs/Rcd > 1-3 can occur in type II 
SNRs with very efficient DSA. The acceleration 
efficiency for the 7/i n j = 10 -4 model in Fig. 10 
(light-weight dashed curve) is greater than 50% 
over most of its 400 yr lifetime. This may be rel- 
evant for remnants like Cas A and 1E0102. 2-7219 
which show Rfs/Rcd ~ 1.4. 

3.3. Radio emission vs. ejecta profile and 
age 

It is well known that young SNRs with power- 
law ejecta and power-law ambient medium den- 
sity profiles have self-similar solutions if CR pro- 
duction is absent or unimportant (i.e., Chevalier 
1982a, b). This will be true for the efficient pro- 
duction of CRs as well if the CR production is 
time invariant (Chevalier 1983). If nonlinear DSA 
occurs and the acceleration efficiency varies with 
time, the self-similarity is broken (see Ellison et al. 
2004), as is the case with an exponential ejecta 
density distribution (e.g., Dwarkadas 2000), or for 
a power-law ejecta distribution once the reverse 
shock enters the plateau region of the ejecta. 

In Fig. 11 we show radio emission profiles at 
various <snr for type la models with 7?i n j = 10~ 3 
having exponential (top panel) and power law 
(bottom panel) ejecta density profiles. In self- 
similar evolution, the ratio Rps / Rcd remains con- 
stant and this is approximately the case for a 
power-law ejecta density profile for isNR £ 300 
yr. At later times, the self-similarity is broken, 
as is the case at all times for exponential ejecta 
density profiles. The light-weight solid curves are 
test-particle profiles at 150 yr for comparison. 

Besides Rfs/Rcd, the structure of the radio 
emission in the interaction region between the CD 
and the FS depends on the assumed ejecta distri- 
bution and on the age of the SNR. At early times 
for the power-law case (solid curve, bottom panel 
of Fig. 11), the radio emission peaks near the con- 
tact discontinuity. This result is consistent with 
the self-similar model described in Cassam-Chena'i 
et al. (2005) but, as discussed above, depends 
somewhat on the starting conditions of the CR- 
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hydro model. At later times the emission drops 
inside the FS and, as expected, the details of the 
ejecta profile cease to be important. The curves 
for 1-10 keV X-rays are not shown, but due to ra- 
diative losses and contrary to the radio, they peak 
strongly just behind the FS for all isNR as shown 
in Fig. 8. 

3.4. Acceleration efficiency 

In Fig. 12 we show the acceleration efficiency, 
i.e., the fraction of energy flux crossing the shock 
that goes into relativistic ions (see eq. 13 of Ellison 
et al. 2000), for various 7yi n j (light-weight curves) 
and the fraction of total SN explosion energy put 
into CRs, Ecr/Esn, forryinj = 10~ 4 (heavy-weight 
dashed curves). These models use our type la and 
II prototype parameters. For the, perhaps, ex- 
treme case of ?7i n j = 10~ 3 , the fraction of bulk 
flow energy flux (in the shock rest frame) that is 
placed in relativistic ions is > 80% during the 1000 
yr span shown for both SNR prototypes. Even for 
?7inj = 10~ 4 , the efficiency is > 10% most of the 
time and more than 10% of the total SN explo- 
sion energy can be put into CRs over the 1000 yr 
lifetime. 

Of course the actual injection efficiency of SNR 
shocks is uncertain and, as noted by Volk et al. 
(2003), injection may vary over the surface of the 
SNR and be significantly less where the magnetic 
field is highly oblique (see Rothenflug et al. 2004, 
for a discussion of parallel versus oblique shock 
geometry in SN 1006). V61k et al. (2003) esti- 
mate that to supply the galactic CR population 
the overall efficiency need only be ~ 20% of the 
maximum values obtained by DSA. Dorfi (1990) 
and Bcrczhko et al. (1996) obtained similar values. 
Nevertheless, if the shocks in supernova remnants 
accelerate cosmic ray ions this efficiently via dif- 
fusive shock acceleration, clear signatures of the 
acceleration will be present in the radiation pro- 
duced by electrons. 

4. DISCUSSION 

4.1. Narrow interaction region 

Perhaps the most unambiguous indication of 
efficient CR production in SNRs is an inter- 
action region between the contact discontinuity 
and the forward shock which is considerably nar- 
rower than predicted without efficient acceleration 



(e.g., Blondin & Ellison 2001). While the ratio 
Rfs I Rcd depends on various parameters, efficient 
DSA can easily result in the FS being less than 
half the distance ahead of the CD predicted with 
test-particle acceleration (see Figs. 3, 4, and 9). 
This may explain observations of -Rfs/-Rcd which 
are considerably less than the smallest value pre- 
dicted by test-particle, self-similar models, as is 
the case for Tycho's (e.g., Reynoso et al. 1997; 
Decourchelle et al. 2001) and Kepler's (e.g., De- 
Laney et al. 2002; Cassam-Chena'i et al. 2004) 
SNRs. 

Even in SNRs such as Cas A and 1E0102. 2-7219 
in the Small Magellanic Cloud (e.g., Gotthclf et al. 
2001; Gaetz et al. 2000; Hughes et al. 2000), where 
the FS and CD are well separated, DSA may be 
quite efficient. As shown in Figs. 9 and 10, moder- 
ately efficient acceleration and/or the presence of a 
pre-SN wind can result in Rfs/Rcb ^ 1-3. Thus, 
while the observation of -Rfs/-Rcd = 1.0 — 1.1 can 
be explained naturally if CR ions are being pro- 
duced efficiently in type la SNe, larger values of 
Rfs/Rcd do not necessarily exclude efficient ac- 
celeration but may be representative of type II SNe 
with pre-SN winds. 

4.2. Precursor and small-scale structure 

In some SNRs extremely small spatial scales 
in X-ray emission are observed at the FS. Us- 
ing Chandra observations, Long et al. (2003) and 
Bamba et al. (2003) have independently examined 
emission profiles in several thin filaments in pro- 
jection in the northeast shell of SN 1006 which 
show scale lengths as short as 0.04 pc (assuming 
a distance to the SNR of ~ 2kpc). 

In Fig. 13 we compare our type la prototype 
model with rn n j = 10~ 3 to the SN 1006 observa- 
tions. We represent the observations with dashed 
lines which roughly span the maximum and min- 
imum scale heights determined by Bamba et al. 
(2003) (see their Table 4). Even though we have 
not attempted a detailed fit to SN 1006, it's 
clear that our compressed B model (solid curve) 
matches the overall observations quite well and 
the shortest scale heights are extremely well mod- 
eled. As emphasized by Berezhko et al. (2003), 
the shortest scale heights occur inside the forward 
shock and are produced by projection effects when 
B is compressed and there is a sharp drop in emis- 
sivity behind the shock. The actual upstream pre- 
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cursor (indicated with a "P" in Fig. 13) has a much 
longer scale height as expected from TeV electrons 
but is not easily discernable with Chandra against 
background emission. 

While our efficient acceleration model with 
compressed B fits quite well, our uncompressed 
model (dotted curve) clearly does not fit, nor does 
a test-particle model (not shown) , as is clear from 
examining the bottom panel of Fig. 7. As far as 
we can tell, our results are in complete agreement 
with those of Berezhko et al. (2003) (see also Ballet 
2005) and provide convincing evidence for highly 
compressed magnetic fields and efficient DSA. 

4.3. Adiabatic and synchrotron losses and 
the offset of radio and X-ray peaks 

Nonthermal X-ray emission in a fixed energy 
band is very sensitive to both adiabatic and ra- 
diative losses. For typical SNR parameters, syn- 
chrotron X-rays are produced in large part by 
the exponential tail of the electron distribution. 
Therefore, any energy loss results in a large drop 
in emissivity. This contrasts with the adiabatic 
losses of the electrons producing radio emission. 
Since radio is produced by lower energy electrons 
in the power law portion of the distribution rather 
than the exponential part, emission in a fixed en- 
ergy band is less sensitive to adiabatic losses. If 
nonlinear effects from efficient DSA arc important, 
the fixed band radio is even less effected by adia- 
batic losses since the portion of the electron dis- 
tribution producing radio is likely to be concave, 
i.e., flattering with increasing energy. 

The synchrotron loss rate will be greater if the 
magnetic field is compressed at the shock and, 
therefore, will depend on the acceleration effi- 
ciency. As we show in Cassam-Chenai' et al. (2005) 
and in Fig. 9 here, the morphology of the X- 
ray emission near the FS varies noticeably with 
77inj , peaking more strongly as the acceleration ef- 
ficiency increases since electrons lose energy before 
convecting far downstream. This feature provides 
an important diagnostic for acceleration efficiency. 

A direct consequence of X-ray emitting elec- 
trons suffering more losses than radio emitting 
ones, is an offset in the peak emission of the pro- 
jected flux at the FS. As shown in Fig. 8, the ra- 
dio emission peaks well within the X-ray emission. 
The separation will depend on the acceleration ef- 



ficiency since, for a given set of supernova param- 
eters, models with efficient DSA have larger com- 
pression ratios and larger downstream magnetic 
fields. The larger the field, the sharper is the drop 
in X-ray emission behind the shock, and the closer 
to the FS position with be the peak X-ray emis- 
sion. 

5. CONCLUSIONS 

We have presented a detailed discussion of the 
influence of efficient diffusive shock acceleration on 
the radial profiles of synchrotron emission in young 
SNRs. The evidence that collisionless shocks, 
in general, can accelerate particles with high ef- 
ficiency is convincing. There are direct space- 
craft observations confirming it (e.g., Eichler 1981; 
Blandford & Eichler 1987; Ellison et al. 1990; Bar- 
ing et al. 1997; Shimada et al. 1999), plasma simu- 
lations show efficient acceleration consistent with 
spacecraft observations (e.g., Scholer et al. 1992; 
Ellison et al. 1993; Giacalone et al. 1997), Galactic 
cosmic-ray energetics and composition suggest it 
(e.g., Axford 1982; Ellison et al. 1997), and theo- 
retical models certainly allow it (e.g., Axford et al. 
1977; Drury 1983; Ellison & Eichler 1984; Jones 
& Ellison 1991; Berezhko et al. 1996; Malkov & 
Drury 2001; Kang et al. 2002; Blasi 2002). An 
unresolved question, of course, is whether or not 
shock acceleration is efficient in SNRs. If DSA 
is as efficient in accelerating ions as suggested, 
the acceleration process will be nonlinear and will 
noticeably modify the SNR structure and evolu- 
tion. We have shown for typical type la and type 
II SN parameters that these structural changes, 
most important of which is the increased shock 
compression, produce clear signatures in the syn- 
chrotron radiation emitted by electrons. We note, 
incidentally, that signatures in the thermal emis- 
sion may also be present since the energy which 
goes into relativistic ions comes out of the bulk 
thermal plasma and produces a drastic reduction 
in the temperature of the shocked gas (e.g., De- 
courchelle et al. 2000; Hughes et al. 2000; Ellison 
et al. 2004). 

Of course, our assertion that the nonlinear ef- 
fects seen in the structure of SNRs are evidence 
for the efficient acceleration of ions rather than 
electrons depends on how the energy of shock ac- 
celerated particles is distributed between electrons 
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and ions. While no definitive theory exists describ- 
ing this partition, the source of the energy going 
into superthermal particles is the bulk kinetic en- 
ergy of the converging upstream and downstream 
plasmas. Diffusive shock acceleration occurs, at 
its most basic level, when particles diffuse across 
the shock and scatter nearly elastically off the con- 
verging plasmas on either side of the shock. When 
particles are accelerated from the thermal back- 
ground, this process favors heavy particles and it is 
generally assumed that shocks put far more energy 
into ions than electrons. There is direct evidence 
for this disparity in acceleration efficiency at the 
low Mach number shocks which have been stud- 
ied in the hcliosphere (e.g., Fcldman 1985; Shi- 
mada ct al. 1999) (see also Ellison et al. 1994), but 
there is no direct evidence, one way or the other, 
in the much stronger shocks which exist outside 
of the heliosphcrc. Nevertheless, with some confi- 
dence, we believe the structural changes we have 
discussed are produced by ion acceleration with 
the radiating electrons being passive markers of 
the effect. 9 

While direct evidence for the production of 
CR ions in SNRs would be the observation of a 
pion-decay spectral feature in GeV-TeV 7-rays, 
such 7-rays are difficult to detect with the signifi- 
cance necessary to distinguish a pion-decay feature 
from inverse- Compton or bremsstrahlung radia- 
tion. Furthermore, in low density regions, inverse- 
Compton may outshine pion-decay emission, leav- 
ing the question of CR ion production for these 
SNRs open regardless of the sensitivity of 7-ray 
telescopes. The best chance of seeing a strong 
pion-decay signal is when a SNR interacts with a 
dense medium such as the synchrotron-dominated 
SNR RX J1713. 7-3946 (G347.3-0.5) interacting 
with molecular clouds (see Cassam-Chena'i et al. 
2004, and references therein). HESS (High Energy 
Stereoscopic System) has recently measured, with 
high significance, the 1-10 TeV energy spectrum 

9 We note that so-called shock surfing has been suggested 
by a number of workers as an effective way of transferring 
shock energy into electrons (see Hoshino & Shimada 2002, 
for example, and references therein). A thorough discussion 
of this mechanism is beyond the scope of this paper, but 
we note that while some descriptions of this effect show 
large energy gains by electrons, nonlinear effects are almost 
certain to limit the effectiveness of this process (see Scholer 
ct al. 2003), particularly in the strong shocks we envision 
for young SNRs. 



in this remnant (Aharonian et al. 2004) and in 
SNR RX J0852.0-4622 (Aharonian et al. 2005) and 
while pion-decay is certainly the most likely emis- 
sion mechanism, it is not possible, based on TeV 
emission alone, to reliably determine the differ- 
ent 7-ray components in these spectra. It should 
now be possible to test for pion-decay emission us- 
ing the morphology since HESS has, for the first 
time, produced 7-ray images of these remnants, 
and the morphology of inverse-Compton and pion- 
decay should be quite different. Observations 
in the MeV range by GLAST should help signifi- 
cantly to distinguish pion-decay from lepton emis- 
sion and may provide incontrovertible evidence for 
or against SNRs as the source of CRs ions. 

We have emphasized here that another signa- 
ture of efficient cosmic-ray ion production is the 
large reduction in the ratio of the radius of the 
forward shock to the radius of the contact dis- 
continuity, Rfs/Rcd- If a large fraction of the 
shock energy goes into relativistic particles and 
high-energy particles that escape from the shock 
system, r to t ^> 4 and the interaction region be- 
tween the CD and FS will be denser and Rfs / Rcd 
will be smaller than with inefficient acceleration 
(Figs. 3, 4, and 9). This effect may explain ob- 
servations of Rfs/Rcd ~ 1 in Tycho's and Ke- 
pler's SNRs. Type II SNe with pre-SN winds 
may experience efficient DSA yet still show large 
Rfs/Rcd ~ 1.3-1.4, consistent with observations 
of Cas A and 1E0102. 2-7219 (Figs. 9 and 10). 
While complicating factors such as an irregular 
ambient medium, dense knots, thin sheets of emis- 
sion, etc., exist in all SNRs, efficient DSA offers 
a natural explanation for this important aspect of 
SNR morphology. Just as important, a large value 
of Rfs/Rcd observed in a type la SNR is strong 
evidence against efficient DSA. 

Yet another sign of efficient DSA is the pres- 
ence of short scale heights seen in nonthermal X- 
ray emission. Short scale heights are predicted 
with efficient DSA because the shock will strongly 
compress the downstream magnetic field and syn- 
chrotron losses will lower the emissivity immedi- 
ately behind the FS. This results in several related 
morphological effects. First, thin sheets of X-ray 
emission (e.g., Fig. 9) should be common at the 
FS, as is consistent with observations. Second, 
projection effects should result in the distinct sep- 
aration of the radio and X-ray peaks (e.g., Fig. 8), 
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also commonly observed. Finally, as we show in 
Fig. 13, the short scale heights seen in SN 1006 
(e.g., Bamba et al. 2003), are most naturally ex- 
plained as sharply peaked emission behind the FS 
seen in projection (Bcrczhko et al. 2003, have al- 
ready concluded this for SN 1006). The actual 
upstream precursor has a long scale length, as ex- 
pected for TeV electrons, but is weak enough to 
avoid detection. 

Supernova remnant SN 1006 seems a clear case 
where the efficient production of CR ions is tak- 
ing place, but remnants such as Tycho's and Ke- 
pler's, with Rfs/Rcd ~ 1) are also likely candi- 
dates. The presence of a significant population of 
CR ions in young SNRs produces effects that are 
readily observable in radiation produced by elec- 
trons and we have made predictions, capable of 
being tested with Chandra and XMM-Newton, to 
test this assertion. 
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Log 10 Momentum [m p c] 

Fig. 1. — Phase space electron distributions mul- 
tiplied by p 4 . In all panels, the dashed curves were 
produced, without losses, in shells at the FS at the 
times indicated, i s hock- The solid curves are these 
distributions at isNR = 10 3 yr with adiabatic and 
synchrotron losses taken into account and, in the 
top two panels, the dot-dashed curves are these 
distributions at £snr = 10 3 yr with only adiabatic 
losses included. The heavy (light) dotted curve 
in the bottom panel is the total electron (proton) 
distribution, with all losses, from the interaction 
region between the CD and the FS at £snr- The 
magnetic field is compressed at the shock as in 
eq. (4) and 77 inj = 10~ 3 . 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.2. 
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Fig. 2. — Precursor flow speed u and magnetic 
field, B, as functions of z, the distance upstream 
from the forward shock. The magnetic field is in 
units of Bo, the far upstream value, and u(z) is in 
units of the forward shock speed, V^ s k, measured in 
the rest frame of the SN. For the type la models 
shown here, we take B = 10 ^G. For our type 
II wind models (not shown), we assume that the 
length scale of the wind is large compared to the 
precursor scale so that B(z) in cq. (9) is obtained 
using a constant B , where B is the immediate 
upstream field at ^snr- 
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Fig. 3. — Radial synchrotron emission in three 
energy bands for two magnetic field configurations 
and two shock injection efficiencies, rji n j. In all 
panels, the solid curve is radio (1-1.4 GHz), the 
dashed curve is low energy X-rays (0.1-1 keV), 
and the dotted curve is high energy X-rays (1-10 
keV). In the top two panels B is compressed as 
in equation (9), while in the bottom two panels, 
B-2 = Bq. The flux of each band is normalized to 
its value at the FS. 
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Fig. 5. — Radial synchrotron emission for 1-10 
keV X-rays. The dotted curve shows the pro- 
file with both adiabatic and synchrotron losses in- 
cluded, the dashed curve has only adiabatic losses, 
and the solid curve was calculated with no losses. 
Since there are no adiabatic losses in the precur- 
sor, the dashed and solid curves are identical in 
the precursor. 



Fig. 4. — Radial synchrotron emission for three 
energy bands as in Fig. 3 for our type II SNR 
prototype. 
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Fig. 6. — Magnetic field strength in the interaction 
region between the CD and FS for a type la SNR 
(solid curve) and a type II SNR (dashed curve). 
The number above each solid dot is the compres- 
sion ratio that parcel of gas experienced when it 
was shocked. The open circles at the ends of the 
curves indicate the unshocked field at the end of 
the simulation, i.e., at tsNR = 400 yr. 
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Fig. 7. — Line-of-sight projections for the radial 
distributions with compressed B shown in Fig. 3. 
In all panels, the solid curve is radio (1-1.4 GHz), 
the dashed curve is low energy X-rays (0.1-1 keV), 
and the dotted curve is high energy X-rays (1-10 
keV). 
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Fig. 8. — Linc-of-sight projections for the radial 
distributions shown in Fig. 3 normalized to the 
forward shock radius. The magnetic field is com- 
pressed in both panels. As in the previous figures, 
the solid curve is radio (1-1.4 GHz), the dashed 
curve is low energy X-rays (0.1-1 keV), and the 
dotted curve is high energy X-rays (1-10 keV). 
Note that the radio emission (solid curves) peaks 
well within the X-ray emission in all cases. The 
fluctuations, most noticeable in the radio emission 
for ?7i n j = 10~ 5 , are numerical noise. 
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Fig. 9. — Light-of-sight projections for the 1-10 
keV X-ray band for various injection efficiencies 
as marked. All results include magnetic field com- 
pression and are for £snr = 400 yr. While the 
absolute normalization is arbitrary, the relative 
normalization between the various plots is correct 
(note the different vertical scales in the two pan- 
els). 
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Fig. 10. — Line-of-sight projections for 1-10 keV 
X-rays for our type II prototype with compressed 
magnetic field. The light-weight dashed curve in 
the top panel has 77i n j = 10~ 4 , while all other 
models shown have 7yi n j = 10~ 3 . In the top 
panel, v w = 20 km s^ 1 and dM/dt (in units of 
M Q yr" 1 ) is varied as indicated. In the bottom 
panel, dM/dt — 10~ 5 M Q yr -1 and v w is varied. 
The light-weight dashed curve in the top panel 
with 77 in j = 10~ 4 has an acceleration efficiency at 
isNR = 400 yr of more than 70% (similar to that 
shown in Fig. 12) and demonstrates that values 
of Rfs/Rcb ~ 1-4, as observed for Cas A, are 
consistent with efficient DSA in type II SNe. As 
in Fig. 9, the absolute normalization is arbitrary 
but the relative normalization between the various 
plots is correct. 
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Fig. 11. — Radio synchrotron profiles for our type 
la prototype at various ages, as indicated, and 
for an exponential ejecta distribution (top panel) 
and a power law ejecta distribution with n = 9 
(bottom panel). The line styles indicate the same 
ages in both panels and in all cases, except for 
the curves marked TP:150yr, fy in j = 10 -3 and the 
ambient magnetic field is B = 10 ^G. As shown 
by the solid curve in the bottom panel, a power 
law ejecta distribution produces a radio profile at 
early times that peaks near the CD. The light- 
weight solid curves are test-particle results shown 
for comparison. 
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Fig. 12. — The light-weight curves are the shock 
acceleration efficiencies for our two SN prototypes 
for various ?7i n j, as indicated. The heavy-weight 
dashed curves in both panels are the fractions of 
total SN explosion energy going into CRs for the 
case where 77; n j — 10~ 4 . The line styles indicate 
the same values of ?7i n j in both panels. 
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Fig. 13. — Comparison of X-ray line-of-sight pro- 
files from the CR-Hydro model with Chandra ob- 
servations of SN 1006. The dashed curves roughly 
span the maximum and minimum scale heights de- 
termined by Bamba et al. (2003) where they as- 
sumed exponential profiles. Using their Table 4, 
we set the maximum (minimum) upstream scale 
height to be 3 (1) arcsec, and the downstream 
maximum (minimum) scale height to be 30 (10) 
arcsec (the radius of SN 1006 is about 0.25°). The 
solid curve is the X-ray emission in the 1.2-2 keV 
band using our compressed B model and for com- 
parison, we show (dotted curve) the 1.2-2 keV 
band without compressing the field. We have po- 
sitioned the peaks of the dashed curves to match 
the solid curve. 
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